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Abstract 

The existence of non-singular string cosmologies is established in a class of two-dimensional 
supersymmetric Hybrid models at finite temperature. The left-moving sector of the Hybrid 
models gives rise to 16 real (A/4 = 4) spacetime supercharges as in the usual superstring 
models. The right-moving sector is non-supersymmetric at the massless level, but is charac- 
terized by MSDS symmetry, which ensures boson/fermion degeneracy of the right-moving 
massive levels. Finite temperature configurations, which are free of Hagedorn instabilities, 
are constructed in the presence of non-trivial "gravito-magnetic" fluxes. These fluxes in- 
ject non-trivial winding charge into the thermal vacuum and restore the thermal T-duality 
symmetry associated with the Euclidean time circle. Thanks to the unbroken right-moving 
MSDS symmetry, the one-loop string partition function is exactly calculable beyond any 
a'-approximation. At the self-dual point new massless thermal states appear, sourcing lo- 
calized spacelike branes, which can be used to connect a contracting thermal Universe to an 
expanding one. The resulting bouncing cosmology is free of any curvature singularities and 
the string coupling remains perturbative throughout the cosmological evolution. 
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1 Introduction 



One of the most acute open problems of Modern Cosmology is that of the initial singularity. 
Assuming that the matter content of the Universe satisfies the weak or the null energy 
condition, and extrapolating the cosmological evolution arbitrarily back in time by using the 
equations of quantum field theory and Einstein's theory of general relativity, one is driven 
to a curvature singularity, the "Big Bang", signifying a break-down of our description of 
the underlying physics. Moreover, the singularity theorems [lj strongly suggest that it is 
unlikely to overcome this problem within the framework of quantum field theory and general 
relativity. 

In string theory, however, new degrees of freedom arise and various non-geometrical 
phenomena take place, opening up new possibilities to address the initial singularity problem 
[2 10 . Indeed, the full set of stringy degrees of freedom cannot be incorporated within a 
single effective field theory description. Rather, there exist distinct effective descriptions valid 
in specific restricted regions of the moduli space. Typically in such regions, most stringy 
states carry large masses of the order of the string scale or higher. At low enough energies, 
these can be integrated out giving rise to an effective field theory in terms of a finite number 
of light fields. However, as the moduli fields vary adiabatically, we encounter points where 
new states become relevant, forcing us to switch from the initial description into another, 
using the duality symmetries of string theory. An apparently singular geometry may then be 
mapped into a non-singular, smooth geometry in the dual description. Furthermore, around 
the transition points a conventional field theory description may be absent, as the full string 
theory dynamics becomes relevant. For comprehensive discussions, see e.g. (ij. 

In a cosmological context, where a number of moduli fields and the background metric 
acquire time dependence, we expect the various dual descriptions to be connected dynam- 
ically via non-trivial phase transitions |2||H|8|[9| [Ilj - [l3] . The cosmological evolution will be 
marked by transition eras, where new stringy states become light and dominate the dynam- 
ics. In general, the transitions among the various dual phases take place at finite values of 
the moduli, where the underlying target space background has finite size and/or curvature. 
Thus, one expects these transitions to occur before the appearance of cosmological singular- 
ities. The hope is that, once the underlying string dynamics is properly taken into account, 
the various non-singular cosmological phases would be connected consistently, successfully 
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covering the whole history of the cosmological evolution. Previous work on string dualities 
and cosmology includes (2 10. 



The greatest obstacle towards a concrete realization of this attractive cosmological sce- 
nario is the appearance of Hagedorn instabilities in string theory at finite temperature and 
in the presence of other supersymmetry-breaking sources. Due to the exponential growth in 
the number of single-particle states as a function of mass, the partition function of super- 
strings at finite temperature diverges once the temperature exceeds the Hagedorn temper- 
ature (Tjj ~ M s ). This divergence, however, is not a pathology of string theory. Rather, 



it signals a phase transition towards a new thermal vacuum 11 13 15 . In the Euclidean 
description of the thermal system, where the time direction is compactified on a circle of 
radius R , T = l/(2nR ), certain stringy states winding the Euclidean time circle become 
tachyonic precisely when the thermal modulus Rq exceeds its Hagedorn value. Thus, the 
Hagedorn divergence can be interpreted as an IR-instability of the underlying Euclidean 
background and the phase transition is driven by tachyon condensation. The crucial point 
is that, in the presence of such condensates, the thermal vacuum acquires non-trivial wind- 



ing charge 11, 13 . The question that we would like to address is whether it is possible to 



identify stable thermal vacua, characterized by non-zero winding charge, whose dynamics 
could be quantitatively treated by utilizing e.g. the powerful 2d worldsheet CFT techniques 
of (weakly coupled) string theory. 

Considerable progress towards this direction has been achieved in the context of (4, 0)- 
supersymmetric type II vacua in arbitrary dimension, where it was found that the intro- 
duction into the thermal system of certain discrete "gravito-magnetic" fluxes, associated 
with the graviphoton (Giq) and the axial vector (Bio) backgrounds, lifts the tachyonic in- 



stabilities 16 17 . Essentially, these fluxes inject non-trivial winding charge into the thermal 
vacuum, and restore the stringy T-duality symmetry R — > R 2 C / Rq of the finite-temperature 
system^] The self-dual point occurs at the fermionic point R c = 1/ \^2, where new massless 
thermal states appear. These states correspond to vortices carrying non-trivial momentum 
and winding charges along the Euclidean time circle, and lead to an enhancement of the 
local U(1)l gauge symmetry to a non-Abelian [SU(2)L]k=2 gauge symmetry. The one-loop 
string partition function is finite for all values of the thermal modulus Rq and reduces to the 



1 In the absence of the gravito-magnetic fluxes the type II thermal partition function does not enjoy 
thermal duality symmetry. In the heterotic case, however, this duality symmetry is present. 



2 



conventional thermal partition function in the large Ro limit. 

In this paper we will study the induced cosmological evolution associated with such ther- 
mal configurations in a class of two dimensional (4, 0)-vacua, the so called Hybrid vacua, 
introduced in [17]. These are very special vacua: the right-moving supersymmetries are 
broken spontaneously at the string level, but are replaced by the recently-discovered Mas- 
sive (boson/fermion) Spectrum Degeneracy Symmetries (MSDS) of [g)[l8]. The presence of 
particular discrete "gravito-magnetic" fluxes renders the finite temperature configurations 
free of tachyonic/Hagedorn divergences and, in addition, the thermal vacuum is character- 
ized by unbroken right-moving MSDS symmetry. As a result of the latter symmetry, the 
one-loop string partition function can be computed exactly in a', with the contributions of 
the massive bosonic and fermionic towers of string states^] exactly cancelling each other. As 
the value of the thermal modulus Ro is varied, the thermal system exhibits an interesting 
phase transition, occurring precisely at the extended symmetry point Rq = l/\/2. At this 
point additional massless states appear, driving the phase transition. The two phases are 
related by thermal duality symmetry Ro — > l/(2i?o) (with the extended symmetry point 
being the self-dual point )|^] In both phases the temperature, the energy density and pressure 
are bounded. They achieve their maximal critical values at the extended symmetry point. 
Thanks to the unbroken MSDS symmetry, the equation of state characterizing each phase 
is effectively that of massless thermal radiation in two dimensions. 

The resulting cosmological evolution is non-singular, and it is marked by a phase tran- 
sition between the two dual thermal phases. We will show that the phase transition can be 
described in terms of a space-like brane, corresponding to a localized pressure that is sourced 
by the additional massless thermal modes that become relevant at the critical temperature. 
The role of the brane is to glue the two dual thermal phases in a consistent way. As we 
will see, both the dilaton and the scale factor bounce across the brane. For the scale factor, 
the bounce corresponds to a smooth reversal of contraction to expansion, with the Big Bang 
curvature singularity of general relativity being absent. There is a discontinuity in the first 
time derivative of the dilaton, determined by the tension of the brane. This discontinuity is 
naturally resolved by the presence of extra massless thermal states localized at the brane. 



2 Notice that, at the fermionic point, one may treat lattice states as fermionic oscillator states in a uniform 
manner. 

3 Various aspects of thermal duality symmetry have been also discussed in 19-21 
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The system remains weakly coupled throughout the cosmological evolution, with the dilaton 
attaining its maximal value at the phase transition. This maximal value is determined in 
terms of the brane tension, the critical temperature and the number of the thermally excited 
massless degrees of freedom in the system. 

We argue that the gluing mechanism provided by such spacelike branes is a fundamental 
dynamical phenomenon in string theory, capable of connecting consistently the space of 
"momenta" with the dual space of "windings". Such branes will typically appear at the 
self-dual points, which are characterized by enhanced gauge symmetry. The presence of 
additional marginal operators with non-trivial momentum and winding charges is crucial in 
order for a dynamical transition between the two dual spaces to occur. 

We would like to stress here that our stringy, singularity-free cosmological solution differs 
from various incarnations of the pre- Big Bang scenario p], which require the application 
of a non-perturbative S-duality symmetry in order to resolve the very early, pre- Big Bang 
phase. In fact, the occurrence of strong coupling dynamics is often another great obstacle in 
realizing the most general cosmological scenario outlined above. It is an open and interesting 
question to find branes that can dynamically glue cosmological phases that are related by 
S-duality symmetry. It is also interesting that the two-dimensional cosmological Hybrid 
models are amenable to perturbative treatment throughout their entire history. 

The plan of this paper is as follows: 

In section [2] we summarize the main features concerning tachyon-free, thermal configu- 
rations of (4, 0)-vacua at finite temperature. In section [3] we focus on the two dimensional 
Hybrid models and analyze in detail the resulting phase structure of the thermal system. We 
proceed in section [4] to study the backreaction of the thermal Hybrid state on the initially flat 
metric and dilaton background and exhibit the resulting non-singular cosmology. In section 
|5j we present more general non-singular cosmological solutions associated to distributions of 
spacelike branes. We end up with our conclusions and directions for future research. 



2 Thermal (4,0)-type II string vacua 

Before describing the two dimensional Hybrid models, we review some general results con- 



cerning thermal states of (4, 0)-type II models 16,17 . Such models correspond to freely 
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acting asymmetric orbifold compactifications, which preserve the 16 real supersymmetries 
(A/4 = 4) arising from the left-moving sector of the worldsheet. The right-moving super- 
symmetries are broken spontaneously at the string level, via the coupling of an internal 
cycle, e.g. the X 9 direction compactified on a circle, to the right-moving spacetime fermion 
number Fr = a. At special values of the internal compactification moduli, we encounter 
extended symmetry points with a non-Abelian enhancement of the local gauge symmetry 



group 16, 17 . In the following discussion, we take the values of the internal moduli to be 
frozen at such an extended symmetry point and the string coupling to be sufficiently weak. 

The finite temperature partition function can be computed via a Euclidean path integral 
on the S 1 (Rq) x Ai manifold, where the Euclidean time direction is compactified on a 
circle of radius Rq. Here Ai denotes the spatial manifold. At least one spatial direction is 
taken to be very large (or non-compact). The conventional thermal ensemble corresponds 
to a freely acting orbifold, obtained by modding out with (—1) F 5 , where S is an order- 
two shift along the Euclidean time circle. The coupling to the spacetime fermion number 
F = a+d is dictated by the spin-statistics connection. The resulting thermal system develops 
instabilities, the Hagedorn instabilities, when Rq < Rh = \/2, where certain stringy states 



winding the Euclidean time circle become tachyonic 11 12 14 . These instabilities signal the 



onset of a non-trivial phase transition around the Hagedorn temperature, see e.g. 1113, 15 



Generically, one expects the tachyon condensates involved in the transition to lead to large 
backreaction, driving the system outside the perturbative domain. A precise quantitative 
analysis of the dynamics is still lacking. At the very least, we can identify a crucial property 
of the would-be stable high-temperature phase: since certain winding modes acquire non 
trivial expectation values, (O n ) 7^ 0, the thermal vacuum must be characterized by non- 
zero winding charge. 

In this paper, we will consider tachyon-free thermal configurations where, in addition to 
the temperature deformation, we turn on certain discrete "gravito-magnetic" fluxes associ- 
ated with the graviphoton (Gjo) and the axial vector (Bio) gauge fields. These are geometri- 
cal fluxes threading the Euclidean time circle together with the internal cycle responsible for 



the breaking of the right-moving supersymmetries 16,17 . Effectively, they describe gauge 
field "condensates" which have locally vanishing field strength but a non-zero value of the 
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Wilson line around the Euclidean time circle^] They correspond to non-trivial, topological 
vacuum parameters, defining (at least perturbatively) super-selection sectors of the thermal 
system. The presence of such fluxes modifies the thermal frequencies and the contribution to 
the free energy of various states charged under the graviphoton and the axial vector gauge 
fields. As the would-be winding tachyons are charged under these gauge fields, for large 



enough values of the "gravito-magnetic fluxes", the tachyonic instabilities are lifted 16 17 
As an illustrative example, consider the case where the X 9 -cycle factorizes from the 



rest of the internal compact manifold 16 . In the presence of non-zero G90/-R9 = G and 
-B90/-R9 = B background fields and finite temperature, the contribution to the one-loop 
partition function of the (X°, X 9 ) sub-lattice is given by 

/in /«'-9 ^ e "# {R%\™ +Tn \ 2 +Rl\m 9 +Grh +T(n 9 +Gn°)\ 2 } ^ e 2inB(m 9 n -m°n 9 ) 



T 2 

m.n 



X 



^m (a+a)+n°(b+b) ^_^yh 9 a+n 9 b+m 9 n 9 ^ 



The Euclidean-time A°-cycle couples to the total spacetime fermion number F and the 
internal spatial A 9 -cycle couples to the right-moving fermion number Fr. The conventional 
finite-temperature system (suffering from the well-known Hagedorn pathologies) corresponds 
to G = B = 0, where the AT°-lattice coupling to the total fermion number factorizes. The 



deformed tachyon-free configuration arises at the point G = 2B = 1 16 . This point is 
very special since it leads to an asymmetric factorization of the original lattice into two 
parts, one of which couples only to the left-moving /^-symmetry charges, while the second 
one couples only to the right-moving /^-symmetries. To see this, shift rh 9 — > m 9 — m° and 
n 9 — > n 9 — n° to define the frame in which the (2, 2)-lattice factorizes into two (1, l)-lattices, 
asymmetrically coupled to the left- and right-moving spacetime fermion numbers Fl = a, 
F R = a, respectively: 

R e -^|m°+n°r| 2 f^rhPa+nPb+rhPn x e -^\m B +n B r\ 2 ,_^m 9 a+n 9 b+m 9 n 9 ^2) 

v rh°,n° v m 9 ,n 9 

In this representation, the X° circle couples to the left-moving spacetime fermion number 
Fi = a, and so the configuration can be described as a freely-acting asymmetric orbifold 
(— 1) Fl Sq of the initially supersymmetric (4,0) model. 



See e.g. 22 for discussions of such configurations in the context of local gauge theories at finite temper- 



ature and 20 23 in the context of string theory. 
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In 17 , general classes of (4, 0)-vacua were studied, including cases where the X -cycle 
does not factorize from the rest of the compact spatial manifold. It was found that the point 
G = 2B = 1 leads to the unique tachyon-free, thermal configuration associated with any 
such initially supersymmetric (4, 0)-vacuum. Moreover, it was established that there exists 
a frame in which the full compactification lattice admits an asymmetric factorization, where 
the (1, 1) sub-lattice associated to the Euclidean time cycle couples only to F L , while the 
internal, spatial sub-lattice couples only to Fr\ 

r m [l'l] = r (1)1) E](i? )®r ((i _ M _ 1) [f](G /J , J B JJ ) . (2.3) 



This result generalizes (2.2). The thermal configurations remain tachyon-free under all pos- 
sible deformations of the dynamical, transverse moduli associated with the internal spatial 
sub-lattice, and of the thermal modulus R [17]. Indeed when G = 2B = 1, the 0$0 8 sec- 
tor, which appears in the spectrum, necessarily carries non-trivial momentum and winding 



charges, and its lowest mass satisfies 17 



1 



for any deformation of the transverse, internal spatial lattice. As a result the lowest-lying 
thermal states coming from the dangerous OsOs-sector are at least massless (for details 
see [17]). 

Thus, the lowest mass in the OsOs sector is never tachyonic but becomes massless when 
Rq = 1/a/2 and for particular values of the transverse moduli in the internal spatial lattice, 
giving rise to enhanced symmetry points. The one-loop string partition function is finite for 
all values of the thermal modulus Ro, and it enjoys the T-duality symmetry: R — > l/(2i? )- 
The self-dual fermionic point Rq = l/\/2 corresponds to the extended symmetry point. 

Let us focus on a special point in moduli where the local Abelian gauge symmetry 
associated to the X° and X 9 cycles is extended to: 

U(1) 2 L x 17(1)! x > ([SU(2) L ] k=2 x U(1) L ) x (U(1) R x [SU(2) R } k=2 ) x • • • . (2.5) 

The first factor and the second factor in each parenthesis on the r.h.s. correspond to the 
X°, X 9 -cycles respectively. This enhanced gauge symmetry allows us to re- interpret the 
deformation in terms of the "gravito-magnetic" fluxes as condensates of vortices carrying 



7 



non-trivial momentum and winding charges along the Euclidean-time and the X 9 -cycles. 
Indeed, the marginal deformation switching on the "gravito-magnetic" fluxes is: 

R 2 9 {G + 2B) (dX°dX 9 ) + R 2 9 (G - 2B) (dX 9 dX°) . (2.6) 

When G = 2B = 1, the second term vanishes. At the enhanced symmetry point, we are 
free to perform an SU(2)l x SU(2)r rotation^] in order to re-express the current-current 
deformation as: 

dX°dX 9 (^°cosX°) x (t/> 9 cosX 9 ) , (2.7) 

where ip° {4> 9 ) is the worldsheet super-partner of X 9 (X^). In other words, the deforma- 
tion injects into the thermal vacuum non-trivial momentum and winding charges along the 
Euclidean-time and the X 9 - cycles. The contribution to the free energy of the massive 
string states is effectively regulated, resolving the would-be divergences. This observation 
illustrates how the Hagedorn instabilities can be lifted once the Euclidean thermal vacuum 
acquires non-zero winding charge. 

As we will see, the (Hagedorn-free) deformed thermal vacua exhibit interesting phase 
transitions. In the cosmological context, such tachyon-free thermal configurations will be 
considered in full detail in a class of two-dimensional (4, 0)-vacua, the so-called Hybrid 



models, introduced recently in 17 . In this class of models, the modified free energy in each 
phase satisfies all the usual thermodynamical properties as a function of the temperature. 
Namely, as the temperature increases, the free energy decreases monotonically up to the 
phase transition and the specific heat is positive. At the phase transition, the temperature 
achieves a maximal critical value. As we will see, the resulting phase structure is not only 
free of the Hagedorn instabilities, but also the induced cosmology turns out to be free of the 
initial "Big Bang" curvature singularity. 

3 The tachyon-free thermal Hybrid model 

In this section we study thermal configurations of the two-dimensional Hybrid models in- 



troduced in 17 . The Hybrid models describe very special (4, 0)-supersymmetric vacua: the 



5 The three [5'C/(2)i]fc = 2-currents are dX°, -0°cosX° and ?/> sinX£. At the enhanced symmetry point 
one may choose any one of the three. Notice that the 5'C/(2)-symmetry currents under consideration also 



utilize the fermionic worldsheet super-partner ip° 12 . 



8 



right-moving supersymmetries are broken spontaneously at the string level via the coupling 
of an internal cycle to Fr, and are replaced by the Massive Spectrum Degeneracy Symme- 
tries (MSDS) of (9,18 . As in type II superstrings, we can define two versions of the model, 
Hybrid A and Hybrid B, depending on the relative chirality of the left- and right-moving 
Ramond sectors. In the maximally symmetric case, the partition functions are given by 



Z 



Vo 



(27T) 



d 2 T 



T 



4(Imr) 



E< 

a,b=0,l 



\a+b+[iab @ \b\ 



E 

7,5=0,1 



0[ 



718 



E( 

a,b=0,l 



Q[a]12 
> =12 



(3.1) 

where \l — in Hybrid B (// = 1 in Hybrid A) and V 2 stands for the volume of the two 
large (Euclidean) directions X^X 1 . In terms of SO(n) left- and right-moving characters, 
the partition function takes the form 

1/24 - ^24 J ^*-E s [T) <j ^ 



(3.2) 



(2tt) 2 7^4(Imr) 2 K ' ^ rf ^ 8V ' ' {{V 8 - C 8 ) ' 

where the upper (lower) entry refers to the Hybrid B (Hybrid A) case. Thus, the left-moving 
sectors are described in terms of the 5'0(8)-characters, as in conventional superstring models, 
and in terms of the chiral ^-lattice. The right-moving sectors are described in terms of the 
MSDS S'0(24)-characters. In particular, the right-moving characters satisfy the identity 

V 2i - S 2A = 24, (3.3) 

exhibiting the Massive Spectrum Degeneracy Symmetry of the right-moving sector. 

The massless content of the Hybrid B (Hybrid A) model consists of 24 x 8 bosons and 
24 x 8 fermions, arising from the V$ V 2 a and 5*8 V 2i (C 8 V 2i ) sectors, respectively. There are 
no massless fermions from the right-moving R sector and the RR fields are massive. The 
model can also be described in terms of a freely-acting asymmetric orbifold compactification 
of the type II superstring to two dimensions. The relevant half-shifted (8, 8)-lattice is given 



by 17 



(8,8) lb \ 



r^ 8 x 0[f] 



= E 

rh 1 ,n J ;/,J=2,...,9 



a/ det G LJ JL(G+ J B) /J (m I +rn I )(m J +fn J ) i^a+n^b+r^rfi) 
' — N8 ' 

r 2 f 



(3.4) 

where all internal radii are at the fermionic pointj^J The modular covariant cocycle describes 

6 Note that in this representation the fermionic point does not simply correspond to Ri = 1/^/2, as would 
be the case in simpler compactifications. The reason is the non-trivial coupling to Fr which asymmetrically 



deforms the lattice. The specific values of the 
given in 17 . 



-moduli corresponding to the fermionic/ MSD S point are 
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the coupling of a single internal cycle of the lattice, which we will henceforth refer to as the 
X 9 -cycle, to Fr and leads to the spontaneous breaking of the right-moving supersymmetries. 
At the MSDS point the Gjj and Bjj tensors take very special values. In fact, the MSDS 
structure characterizing the right-moving sector gives rise to an enhanced local non-Abelian 
gauge group: U(l)\ x [SU {2) R ]\ =2 . A detailed study of the classical moduli space of the 
MSDS and Hybrid vacua, as well as the study of marginal deformations connecting these 



vacua to four dimensional supersymmetric vacua, can be found in 17 



At finite temperature, the backreaction on the initially flat metric and dilaton background 
induces a cosmological evolution. As we saw in section [2j thermal tachyon-free configura- 
tions, able to bypass the Hagedorn instabilities, can be constructed in terms of freely-acting 
asymmetric orbifolds along the Euclidean time circle. In particular, these are obtained by 
modding out with (—1) Fl 8o, where Sq is a Z2-shift along the Euclidean time circle, and can 
be interpreted as discrete deformations of the conventional finite temperature models. An 
important result is that they remain tachyon-free under arbitrary marginal deformations of 



the transverse dynamical moduli associated with the compact internal eight-manifold 17 



We will describe below some of their main features. Our notations concerning the wrapping 
numbers (m, n), as well as the the left- and right-moving momenta (pl,Pr) of the r^^-lattice 
are as follows: 



r 2 

m,n m,n 



r„n = ?^?^ 1 g = e 2 ^, and Plr = l=(^±nR^j . (3.5) 
For sufficiently weak string coupling, the thermodynamic behavior of the system is cap- 



tured by the one-loop partition function. In the Hybrid B case, this is given by 16, 17 

Z f d 2 T 



'24 J 



\ - I 8 7r( 'I mr )3/2 ~§T ( V 8 T rn,2n + 8 T m+ i i2n+1 - S 8 T m+ l j2n - C 8 r mi2 „ + l) (V 2i ~ S 2 

(3.6) 

where V\ = 2ttR\, whereas the Hybrid A case is obtained by the chirality exchange Ss —> Cs- 
The (— 1) Fl 8q orbifold acts thermally on the left-movers and breaks the (4, 0)-supersymmetries 
spontaneously. On the other hand, the right-moving MSDS structure remains unbroken. In 
the odd-winding sector, the left-moving GSO projection is reversed, and so the 0%- and C%- 
sectors appear in the spectrum. As the right-moving sector, however, begins at the massless 
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level, the model remains tachyon-free for all values of the thermal modulus Rq (as follows 
by level-matching and the properties of the T m+ i 2n+1 -lattice). The Hybrid B and Hybrid A 
models are mapped into each other under the T-duality transformation R Q — > 1/(2R ). The 
thermal quanta of Hybrid A are mapped into vortices, carrying non-trivial winding number, 
in the Hybrid B model and vice-versa. 

As we already stressed, the one-loop partition function is finite for all values of the thermal 
modulus i?o- In the Hybrid model, the partition function can be computed explicitly thanks 
to the MSDS identity V24 — §24 = 24. The result is exact in a', and it is given in terms of a 



surprisingly simple expression 17 



Rn 



(3.7) 



2R 

The essential feature is a discontinuity in the first derivative of Z as a function of Rq, 
signaling a phase transition. The discontinuity occurs at the self-dual point Rq = 1/ \/2. It 
is sourced by the lowest mass states in the 8 ^24- sec t° r > which become massless precisely at 
the fermionic point. Indeed, the lowest mass in the OsVJm sector is given by 

<v = [Ro ~ , (3-8) 

and it is everywhere positive except for the fermionic point, where it vanishes. Similar 
behavior is exhibited in thermal configurations associated to non-critical heterotic strings in 



2d 24 



At this point, several comments are in order: 

• For Rq > 1/ \/2, the partition function is given by 

Z 24 . . 

v, = K (3 ' 9) 

coinciding with the conventional thermal partition function of 24 x 8 massless bosons and 
24 x 8 massless fermions in two dimension^} The corresponding temperature is T = 1//3, 
where /3 = 2ttR . This result can be obtained by Poisson re-summing over the momentum 
quantum number m, and then mapping the integral over the fundamental domain to an 



integral over the strip 25 , involving the zero winding (n = 0) orbits 8 ! The latter are 



The contribution of each massless boson/fermion pair to the one-loop thermal partition function in 2 
dimensions is 1/(8Rq). 

8 Due to the initial (4, 0)-supersymmetry, the contribution of the modular invariant (m, n) — (0,0) orbit 
vanishes. 
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contained in the even winding sectors VgV^, SsV^, ^s^W and SsSW- The splitting of the 



integral (3.6) into its orbit by orbit contributions, and the subsequent mapping to the strip, is 
only valid for R > l/y/2, since then the individual (m, n) 7^ (0, 0) orbit contributions in the 
Lagrangian lattice are exponentially suppressed as r 2 — >• oo and we can safely exchange the 
orders of summation and integration^} For R = 1/ y/2, the lack of this absolute convergence 



due to the presence of additional massless states, results in the non-analytic part of (|3.7|) 
So for Rq > 1/ \/2 we obtain 



17 



24 



17 



T7 = RoJ2 / o n v e ~^~ ~% ( v * - (-T s 8 ) (V 24 - S 24 ) . (3.10) 
V i J\\ 8vr(Imr) 2 r/ 8 

The sectors VsV^a and SsVza, which contain the initially massless bosons and fermions, are 
thermally excited as in the conventional thermal deformation. The deformation deviates 
from the conventional one in the sectors VgS^ and SgS^zi- These two sectors are massive, 
with masses of the order of the string scale, even at zero temperature. Due to the unbroken 
right-moving MSDS symmetry, V24 — ^24 = 24, and level-matching, the contributions of 
all massive boson and fermion string oscillator states cancel. Only the thermally excited 



massless states give non- vanishing net contributions, leading to the result (3.9): 



1 = (24 x 8) -L- y , 1 X9 = — . (3.11) 
Vi 1 J vr^o ^ (2m + 1) 2 ,Ro 1 } 

• The integral over the strip allows one to rewrite the partition function Z in terms of (the 
logarithm of) a (spacetime) "dressed" thermal trace over the 2 nd quantized Hilbert space 
of the initially supersymmetric Hybrid B model. Because of the asymmetric nature of the 
deformation, this is given by the right-moving fermion index: 

Z = lntr[e-^ (-1) Fr ] • (3.12) 

This expression is strictly valid for (3 = 2nR > 2it/ y/2, as indicated by the non-analyticity 



of the exact result for the partition function (3.7). The right-moving fermion index differs 



from the canonical thermal trace in the odd-F^ sector (a — 1). Such states originate from 
the massive sectors VgS^ and SgS^ (with masses of the order of the string scale). As we 
have seen, in the presence of unbroken right-moving MSDS symmetry, the contributions 



9 For R < 1/V2, we can obtain the result either by applying T-duality or by first Poisson re-summing 
over the winding number n, which effectively maps the system on a large radius. See the discussion below. 
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of all states involving massive string oscillators cancel and so the right-moving fermion 
index reduces to a canonical thermal ensemble over the initially massless fields^} We can 
also understand the cancellations as follows. The contribution to Z from a single bosonic 
(fermionic) field of mass M and right-moving fermion number a is given by 

T ^ln(l^(-l) a e-^), u k = Vk 2 + M 2 . (3.13) 

k 

This shows that the contribution of a bosonic field with right-moving fermion number a 
cancels against the contribution of a fermionic field with the same mass and right-moving 
fermion number 1 — a. The right-moving MSDS symmetry requires that such massive bosons 
and fermions come in pairs. 

• Consider now the canonical ensemble, defined in terms of the conventional (undressed) 
thermal trace. For sufficiently large /3, the contributions of all states involving massive string 
oscillators are suppressed and, effectively, the thermal configuration reduces to a system of 
thermal radiation associated with the initially massless states in the sectors VsV^ and Ss^24- 
In this sense, the deformed model does not differ appreciably from the conventional thermal 
system for large 0. For values of /3 close to the string length, massive string oscillators 
become relevant. The canonical ensemble fails to converge for < (3 H = (2?r)V2, due to 
the exponential growth of the number of single-particle oscillator states. The introduction 
of the discrete gravito-magnetic fluxes effectively regulates this behavior. In their presence, 
the temperature can be further increased up to the fermionic point (5 = (2ir) / v2, until one 
encounters a phase transition, which, in the present model, can be followed perturbatively. 

• For Rq < 1/ v2, the partition function is given by 

§- = 24 x (2Ro). (3.14) 



This result can be obtained by applying the Rq — > 1/(2R ) duality to equation (3.9). Al- 
ternatively, it can be derived by Poisson-resumming over the winding number n in equation 



(3.6) and unfolding the integral over the fundamental domain to an integral over the strip, 
as before. Essentially, T-duality maps the Hybrid B system at small radius Rq to a thermal 
Hybrid A system at a large radius 1/(2R ). The light states occur in the VsVk and CsVz^ 
sectors. In both of them, these are purely winding states (or vortices). They can be inter- 
preted as ordinary thermal excitations, with temperature given by the inverse period of the 
10 It would be interesting to see if these cancellations persist at the higher genus levels. 
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T-dual circle: T — 1/(3, (3 = 2h/{2Rq). This definition of the temperature in the small radius 
regime is in accordance with our expectation from conventional thermodynamics: typically 
the one-loop thermal partition function decreases as the temperature decreases. We see 
that the Hybrid B system at small radii is again effectively cold and better understood as a 
Hybrid A thermal system at large radii. 

• At the self-dual point Rq = l/\/2, there are 24 complex (or 48 real) additional massless 



states arising from the OgV^-sector (see Eq. (3.8)). These states carry simultaneously non- 



trivial momentum and winding quantum numbers along the Euclidean time circle. Notice 
that, at the fermionic point, the corresponding values of the left- and right-moving momenta 
are given by 

p L = ±l, p R = 0, (3.15) 
and there is a further enhancement of the local gauge symmetry group associated with the 
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U(1) L x U{1) R [SU(2) L ] k=2 x U(1) R . (3.16) 



The effect of these states in the partition function (3.7) is to induce the non-analytic part. 
The latter can also be expressed as —2A\m y\, which can, in turn, be understood in the 
following way. Near the extended symmetry point the contribution of the lowest mass states 
in the OgV^-sector is given by 



oo 



48 x / — r 1 / 2 e-^lv 1 = 48 \m oV \ J_ T ^ y -W e -v. ( 3 .17) 



l 4 ?rt 4^ J \ m 2 y 



ov 



In the limit m ov — > 0, the leading m G y-dependent contribution is 

-VP 

40F V 2 



48 \m oV \ -^=Y (-1 ) = -24 \m oV \. (3.18) 



In summary, the system defined in Eq.(3.6) has three characteristic regimes, where dif- 



ferent light states dominate the underlying effective field theory dynamics as the value of 
the thermal modulus is varied. There are distinct Hybrid B and Hybrid A phases related 
by "thermal duality symmetry", R — > 1/(2R ) together with S 8 -H- C 8 , and a phase tran- 
sition between them at the self-dual point R = l/y/2. In the Hybrid B phase R > l/y/2, 
the light thermal states belong to the space of momenta, whereas in the Hybrid A phase 
Rq < l/y/2, the light thermal states belong to the space of windings. The third intermediate 
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regime corresponds to the neighborhood of the fermionic point Rq = 1/ \/2 and governs the 
transition between the two dual phases. 

Thus it is necessary to utilize three distinct local effective field theories in order to 
describe fully the stringy thermal system. Alternatively, one can use a single non-local 
description in which the spacetime variables are doubled, utilizing constraint fields that 



depend simultaneously on variables conjugate to the momentum and winding modes 26 . In 
our specific model, the transition between the momentum and winding spaces can be exactly 
resolved in terms of localized branes, as we will discuss later. 

At the fermionic point, the massless states form a localized [SU (2)i\k=2 symmetric space. 
Operators associated with the extra massless states, which are localized at the fermionic 
point Rq = 1/ \/2, induce transitions between the purely momentum and winding modes. 
To see this, recall that the corresponding left- and right-moving momentum charges are 



O- 



Cs; -1,1 ) ■ (3.19) 



given by equation (3.15). We denote the corresponding vertex operators, integrated over 
the worldsheet super-coordinates, + and 0_ respectively. Next, consider a level-matched, 
purely momentum state in the SsV^ sector. At the critical point, such states with lowest 
mass have Pl = Pr = ±f- Let us focus on the state ||, |). Acting on it with 0_ lowers pi 
by one unit and leaves pr unchanged. The new state has pr = —pi = § and it is a purely 
winding state. It corresponds to a level-matched state in the C 8 V^-sector. Pictorially, the 
action of the 0_ marginal operator can be described as 

1 1 

) -4 ( \ : - . 

2' 2 

The explicit mapping can be understood as follows. In the 0-ghost picture, 0_ corresponds 
to the zero-mode of the holomorphic current J_ = il)°e~ lX ° L . Moreover, the (— |)-picture 
vertex operator associated to the massive spacetime fermion state \S 8 ; |, |) takes the form 
e _< ^/ 2 5'io , a e5 x L+k x< k V 2 4, where S\o ia is the spin-field transforming in the spinorial represen- 
tation of SO (10) (of definite chirality) and V24 is the right-moving oscillator contribution 
from the vectorial representation of SO (24). In our conventions, the spinorial representa- 
tion comes with odd chirality, so that the overall vertex operator be consistent with the 
GSO-projection. The action of J_ on \Ss] |, |) yields: 

L(^V* W W = — i^ 10 ,ei^ + ^^ 24 W + reg., (3.20) 

z — w ap ,M 

where 7 is the familiar Dirac 7-matrix in 10 dimensions. As a result of this action, the 
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chirality of the spinorial representation changes. The new vertex operator then corresponds 
to a physical state in the odd winding sector of the theory. Taking the contour integral picks 
up the simple pole contribution from the OPE and we readily verify that 0_ indeed maps 
the purely momentum state, \S&; |, |), into the purely winding state, |Cg; — |, |), as stated 
above. 

To make the duality properties of the system more transparent, it is convenient to express 
the thermal modulus in terms of the variable a G (—00, +00): 

R = e". (3.21) 

Then the physical, duality-invariant temperature can be written neatly as follows: 

\/2 

T = T C e"H , T c = |- . (3.22) 

This expression makes manifest the fact that the temperature of the system is bounded 
from above, never exceeding a critical value: T < T c . As we vary the thermal modulus 
cr adiabatically from —00 to +00, the system heats up to the critical temperature; it then 
undergoes a stringy phase transition and, subsequently, it cools down. In terms of the 



modulus a and the physical temperature, the partition function (3.7) can be written as 

^ = (24V2) e- |CT| = AT, A = = 48tt, (3.23) 

Vi J-c 

where we introduced the thermal parameter A being proportional to the number of thermally 

excited light states. The two dual thermal phases occur at cr > and a < 0, respectively. 

In each phase the pressure and energy density are given by 



2 



P = AT 
dP 

p = -P + T— = AT 2 . (3.24) 

Therefore, in each phase the equation of state is, effectively, that of thermal massless radia- 
tion in two dimensions: p = P. At low temperature, the two dual phases are distinguished 
by the light thermally excited massless spinors: in the Hybrid B phase these transform under 
the Sg-spinor of the internal S"0(8)-symmetry group and in the Hybrid A phase in terms 
of the conjugate Cs-spinor. In addition, the stringy thermal system is characterized by a 
critical temperature T c and, consequently, by a critical energy density given by: 

p c = A T 2 C = 24/tt. (3.25) 



16 



In the stringy thermal system, new light modes appear at T c , which are responsible for 
the transition between the two dual phases. This transition is implied by the underlying 
T-duality of the model and the existence of an extended symmetry point. In particular, it 
should persist to all orders in the string genus expansion. As we will show in the following 
section, this phase transition admits a brane interpretation that glues together the space 
of "momenta" and the dual space of "windings" . This connection is absolutely transparent 
in the two-dimensional Hybrid model under consideration. In higher dimensional thermal 
models, we also expect similar gluing connections between the spaces of "momenta" and 
"windings", described via localized branes at T c . The simple conical structure of equation 



(3.7) will be replaced by a more involved structure according to the dimensionality of the 



Euclidean spacetime. 

4 Non-singular cosmology 

In this section, we will analyze the backreaction of the tachyon-free thermal Hybrid con- 
figuration on the initially flat metric and dilaton background. In several models the in- 
duced cosmological evolution in the intermediate regime T <C Tjj has been extensively 



analyzed 27 , 28 , where an interesting attractor solution corresponding to a radiation-like 



era was found 28 . Other work in the context of hot string gas cosmology includes 29 . In 
this work, however, the important high temperature regime and the structure of the stringy 
phase transition at the Hagedorn temperature will be analyzed. As we have seen in section 
[3j the thermal stringy system has three characteristic regimes associated to the space of light 
momenta (cr > 0), the space of light windings (a < 0) and the extended symmetry point 
which is localized at a = and connects the two. At sufficiently low temperature, the first 
two regimes are effectively described by local field theories involving the light states of the 
Hybrid B and A models, respectively. At the extended symmetry point (a = 0), the effective 
description is in terms of a one- dimensional non-Abelian gauge theory, based on the gauge 
symmetry group 

H L xH R = (SU(2) x U(l) 8 ) L x (17(1) x SU(2) 8 ) R , (4.1) 

in contrast to the other two regimes where Hl is broken down to the Cartan subgroup 
H L — > (U(l) x U(1) 8 )l. In all regimes the massless states are scalars transforming in the 
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( adjoint, adjoint) representation of Hl x Hr. 

The Hybrid A and Hybrid B phases are mapped into each other by thermal duality, 
implying the presence of a maximal critical temperature, T = T c e~'°"L The non- vanishing 
energy density and pressure induce a cosmological evolution, during which the thermal mod- 
ulus a becomes a monotonic function of time: a — > a(t). Since a is a monotonic function, it 
can be used as a parametric time variable, scanning all three effective field theory regimes 
of the stringy thermal system. The cosmological evolution will be marked by a transition 
from the Hybrid A phase to the Hybrid B phase which, at finite coupling, is well localized 
at cr = 0. The phase transition is driven by the extra thermal states that become massless 
at this point and can be effectively described in terms of space-like branes, localized in time. 



Stringy space-like brane configurations in different context have been discussed in 30 . As 
we will see below, the induced cosmological evolution turns out to be non-singular: thanks 
to the presence of the space-like branes, the "Big Bang" singularity of general relativity is 
absent. 

The scope of this section is to establish the existence of an effective Lorentzian action, 
able to describe the three stringy regimes simultaneously, which goes beyond the leading 
cv'-approximation and is valid up to the genus- 1 level. Our guidelines are the symmetry 
properties and behavior of the thermal Hybrid model and its partition function, namely: 



The exact expression (3.7), beyond any cv'-approximation, of the one-loop partition 
function, exhibiting a specific conical structure that is sourced by the additional mass- 
less states at cr = 0. 

The effective equation of state in both Hybrid A and B phases, which is that of thermal 
massless radiation in 2d. Thanks to the unbroken right-moving MSDS symmetry, the 
contributions to the partition function of all massive string oscillator states cancel. 

The thermal duality relation between the Hybrid A and Hybrid B phases, a — > — cr, 
interchanging simultaneously the left-moving spinors: S$ <H- C&. This symmetry, to- 
gether with CPT invariance (always present in an effective field theory), imply that 
below the critical temperature the corresponding effective actions are related by time 
reversal, namely a(t) — > — c(— t). 

The existence of the extended symmetry point at cr = and its brane interpretation, 
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gluing the two phases A and B. 



The above ingredients lead to an effective two-dimensional dilaton-gravity action, able to 
describe simultaneously and in a consistent way the three regimes of the stringy Hybrid 
model. The effective action is naturally separated into three parte 



n 



S = S + S 1 + <S brane , (4.2) 

where 

(4.3) 



S = J d 2 x e- 2 f ^g(±R + 2(V0) 2 ) 



S 



brane 



Si = I d 2 x P, (4.4) 

— K J dx x da e- 2 * y/g^ 5(a). (4.5) 

- So is the familiar genus-0 dilaton-gravity action (written in the string frame). 

- Si is the genus-1 contribution of the thermal effective potential —P. The absence of a 
dilaton factor in this term indicates its origin as a 1-loop effect. The energy density and 
pressure are given by: 



P = P=~ 2 , /? = /3 c e H , (4.6) 



A 

where f3 is the inverse proper temperature, (3 C = (27r)/v / 2 the inverse critical temperature 



and A = A8n is the thermal parameter computed previously in Eq.(3.23). We emphasize 
once more that our one-loop computation of the energy density and pressure is exact in 
a' as the contributions of the massive string oscillator states cancel due to the underlying 
right-moving MSDS symmetry. 

- <Sbrane * s ^ ne s P ace hke brane contribution at the phase transition. It gives rise to localized 
negative pressure and is sourced by the additional 24 complex massless bosons, Xi = 
1, . . . ,24), at the extended symmetry point a = 0. Indeed, the microscopic origin of this 
term follows from the underlying description of the system at the extended symmetry point, 
giving rise to a localized action: 



5 branc = jdx ^gTie *{-g —^j 



(4.7) 

at cr(x°)=0 



11 In principle, higher-derivative terms have to be included in this action together with the contribution 
coming from the higher genus levels. As we will see later, such terms are suppressed during most of the 
cosmological evolution. In the region close to the phase transition their suppression is controlled by the value 
of the brane tension. The existence of an exact CFT description at the extended symmetry point allows us 
to have quantitative control even in this region. 
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We will be interested in spatially homogeneous solutions for gn and <p. Then the Xi equation 
of motion yields: 

-^ = Xi = OLi + 7iv^n x 1 , (4.8) 



dx 1 

where a, and 7$ are integration constants (invariant under a^-reparametrization). Substi- 
tuting the above solution back into the action, determines the value and sign of the brane 
tension in terms of the "rapidity" coefficients jf 

24 

« = $^N 2 = 24| 7s | 2 . (4.9) 

In the most symmetric configuration all the rapidity coefficients are equal, (7$ = j a , i — 
1, ...,24), and so the brane tension is proportional to the number of the extra massless 
states. 

Generically, in the effective action, several other light fields appear; namely, there are in 
total 64 moduli fields parametrizing the coset manifold 

M = SO(8 ' 8) (410) 
SO(8) x SO(8) y ' 

In all field theory phases of the thermal Hybrid model, the one-loop effective potential lifts 



these flat directions 28 , 31 and the above moduli are attracted to the enhanced MSDS sym- 
metry points. This justifies our inclusion in the action of only the dilaton-gravity background 
fields and the contribution of the brane. Thus, in all subsequent analysis, we will work in 
the case where all remaining dynamical (transverse) moduli are frozen to their MSDS values 
so that the right-moving MSDS symmetry of the theory is preserved. 

We are now in a position to determine the induced stringy cosmological evolution by 
deriving the equations of motion in the presence of the brane term. We parametrize the 
metric as follows (x° = t): 

ds 2 = -N 2 dt 2 + a 2 dx l2 , (4.11) 

where the lapse function N and scale factor a = e x are functions of the time coordinate only. 
Without loss of generality, we choose the phase transition at a = to occur at t — 0. In 
terms of these variables, the action, valid for all times, is given by: 

S = Jdtdxe~ 2 ^ +X (Jy 20A-20 2 - k«J(*)J + JdtdxNe x P. (4.12) 
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Varying with respect to N, 4>, A, we obtain the following system of equations: 

2 - <M = Ve 2 <V, (4.13) 



A - A 20-A + — = N 2 e^P, (4.14) 




- 20 2 + ( A - | ) = ±iW*(P - p) - l -Nn5{t) = -^NkS®. (4.15) 



In deriving the first equation, we use the thermodynamical identity 



12 



dP OP 

following from the fact that in any frame the inverse proper temperature, or the period of 
the Euclidean time direction, is proportional to the lapse function N. The r.h.s. of equation 



(4.15) was simplified further by using the equation of state p = P. The crucial effect of the 
brane is to induce a discontinuity in the first time derivative of the dilaton at t = 0. 

The thermal entropy (per co-moving unit cell) is given by 

S = af3(p + P) = 2A|. (4.17) 

On each side of the brane this quantity is conserved, as can be derived by using the equations 
of motion. At the brane, the temperature is fixed attaining its maximal critical value, /3 = (3 C . 
Continuity in the scale factor across the brane implies the continuity of the entropy and vice- 
versa. Equivalently, the stringy phase transition does not produce latent heat and is similar 
to a second order transition. More precisely, by utilizing the equations of motion, we derive 
the following equation: 

p + X(p + P) = ^e- 2 ^ (A - 20) 6{t), (4.18) 



where the r.h.s. represents the brane contribution. Integrating (4.18), we find that the 
thermal entropy is given by 

S(t) = S(t ) + ^^e- 2m [A(0 + ) + A(O-) - 2 (0(0+) + 0((T))] 0(f), (4.19) 



12 This equation exemplifies a deep connection between gravity and thermodynamics, since in the Euclidean 
description, where N — e' CT ' (with x° ~ x° + 2n), the thermodynamical identity is simply derived by the 
Euclidean gravity equations of motion. This non-trivial connection of gravity and thermodynamics will be 
explored in future work. 
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where $(t) is the step function and to < 0. Continuity of the thermal entropyp] then imposes 
the following condition for the discontinuities in the first time derivatives of the scale factor 
A and the dilaton (f>: 

A+ + A_ = 2 + 0_) , (4.20) 

where from now on we simplify the notation to A(0 =t ) = X±, (f)(0 ) = (f>±- When the thermal 
entropy is conserved, we obtain the following relation between the scale factor and the 
thermal modulus a: 



A = lerl + In 



2A ' 



(4.21) 



Therefore, the scale factor a at any time must be greater than (or equal to) a certain critical 
size 



S 

2A' 



(4.22) 



hinting at the absence of the conventional Big Bang singularity of general relativity. 

The equations of motion simplify considerably by choosing an appropriate gauge. To this 
end, we utilize time-reparametrization invariance and work in the conformal gauge N = a. 
The equations of motion then become: 



(i) : 2(f) 2 -2<f)\ = Ce 2 < t> , 

(ii) : A - 20A = Ce 2 *, 



in 



40 2 = - K e x °5{t), 



(4.23) 
(4.24) 
(4.25) 



where Qq = e 



Ao 



,A(0) 



f3 c S/2K is the value of the scale factor when the system reaches its 



critical temperature and 



C 



4A 



(4.26) 



is a constant being proportional to the number of light degrees of freedom. 

Subtracting Eq.(4.23) from (4.24), we obtain A — 2(f) 2 = 0. Further subtracting the latter 



from Eq.(4.25), we obtain 



A 



ne $(t). Observe that both of these equations 



as well as Eq.(4.25) are independent of the constant C. Their structure clearly shows that 



only the first time derivative of the dilaton is discontinuous across the phase transition, while 
its second derivative gives rise to the delta function brane source. The first time derivative 



13 Some generalized solutions will be discussed in Section [HJ 
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of the scale factor is continuous across the brane: A_ 
in terms of the brane tension: 



A_. The discontinuity in <ft is given 



1 



e Xo K. 



Using condition (4.20) together with the continuity of A, we obtain A(0) 



+ 



conservation of the entropy (4.17) then implies 

A = ^ 



(4.27) 
. The 

(4.28) 



showing that the first time derivative of the temperature is also continuous across the critical 
point. The temperature attains its maximal (critical) value T c at t = and, thus, its 
derivative and the derivative of the scale factor must both vanish at this point: /3(0) = 
A(0) = 0. Eq.(4.20) then implies that the first time derivative of the dilaton flips sign 



across the brane. The net result shows that the dilaton bounces during the 



transition. 



Solving Eq.(4.25) we obtain the following expression for the time dependent dilaton: 

p2</> 



,2#t) 



1 + 



t 



(4.29) 



where 



is a positive constant determined in terms of the brane tension k 

: 1 



e Xo K. 



(4.30) 



The above expression shows that the string coupling <?^ r (i) remains smaller than a maximal 
value, <7g tr (0) = e 2 ^ , for all times and, therefore, the perturbative validity of the model is 
guaranteed, provided that <7s tr (0) is sufficiently small. The dilaton reaches its maximal value 
at the critical point and then decreases monotonically. 

It turns out that the string coupling g s tr(0) is not an arbitrary parameter. Indeed, from 



Eq.(|4.23|) evaluated at the critical point we obtain a relation between C and the derivative 

C 



of the dilaton 



Combined with (4.30), it yields: 



8A 



3vr|7 s 



(4.31) 



(4.32) 
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where the last expression in terms of the rapidity coefficient 7 S is obtained by utilizing 



Eq.(4.9), relating the brane tension k with the rapidity, the value of the critical temperature 



/3 C 1 and the numerical value of the thermal parameter A. The validity of the perturbative 
regime is ensured by the smallness of the rapidity coefficient 7 S . 



The evolution of the scale factor is obtained by solving Eq.(4.24). The result is: 

=2^-1*1 



a\t) 



3 2A 



(4.33) 



l+20_|t| 

It is convenient to perform a coordinate transformation r = 2<p_t, x = 2<j)_x 1 . In the new 
coordinates, the cosmological solution for the metric and dilaton is: 

4 el r l 

K 2 1 + \l 



ds' 



9 s 2 tr( T ) 



( -dr 2 + dx 2 ) 



7TK 



1 



1 + \t\ 192 1 + \r\ 
The metric has no essential singularities; both the Ricci scalar 



(4.34) 
(4.35) 



R 



k 2 e-H 



(4.36) 



4 l + br 

and the string coupling are bounded from above. Their maximal values occur at the critical 
point and are both proportional to the square of the brane tension, 0(n 2 ). As a result, all 
higher-derivative terms, as well as higher-genus corrections, are controllable by a perturbative 
expansion in terms of k, provided that k = 24|7 S | 2 <C 1. 

The above solution, even for large times, differs from the naive radiation-dominated 
cosmology, mainly because of the presence of a running dilaton. To see this, we switch to 
the cosmological frame, where the metric takes the form 

ds 2 = -d£ 2 + a 2 (Odx 2 . (4.37) 



The cosmological time £ is given in terms of the conformal parametric time r as: 
l 



2 T 



,\t\u/2 



du 



T\U 



8tt 




erfi 



1 

71 



(4.38) 



with erfi(x) = |erf(ix) being the imaginary error function and x = x/2. For large cosmo- 
logical times £, the asymptotic properties of the error function imply: 



4 e H/2 

K 



l + O 



:,\r\ 



4 



HO 2 + 



(4.39) 
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In the very late cosmological time regime, the asymptotic expansions for the scale factor and 
string coupling are given by 

/ 1 \ 1 192 

a « ) = l H 1 -tawF + "'J ' i; = rf("«' ! + " lln ( K « 2 + C0 " st ' + '")' 

(4.40) 

For large \k£\ ^> 1, the cosmology asymptotes to the flat thermal Milne Universe. There are 
logarithmic corrections due to the running dilaton, giving rise to the above behavior. For 
small cosmological times |/c£| <C 1, the scale factor and the string coupling have the following 



expansion: 



l + ¥-^«) 2 + ^(l<| 3 ) 



(4.41) 



illustrating very clearly that the cosmological scale factor bounces at the location of the 
brane. The dilaton, on the other hand, exhibits a conical structure, which results in a 
discontinuity in its first time derivative. 



5 Generalized solutions 

There are two linear combinations of the equations of motion which are independent of the 
genus-1 thermal contribution, namely 

20 - 40 2 = -Ke Xo 5(t) , 20 - 2A = -Ke A " 5(t) . (5.1) 

The third equation 

20 2 - 20A = Ce 2 ^, 

depends on the thermal contribution through the conserved quantity C, which is deter- 



mined by the thermal entropy; see Eq.(4.26). The most general solution of the first two 
C-independent equations is the following: 

2<f> p A\t\+Yt 

^ = ?n TT > e2A = e2A ° TT with A = — e x ° . (5.2) 

1 + A\t\ +Bt ' 1 + A\t\ +Bt 2 v ; 

It is very important that the integration constant A, multiplying the conical part \t\, is 
fixed by the brane tension k and is, therefore, positive. All solutions satisfying \B\ < 
A avoid the infinite coupling constant singularity and string perturbation theory remains 
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valid throughout the cosmological evolution. Substituting the general solution to the third 
equation, yields the following relations among the integration constants: 

(A + B)(A + T) = 2C + e 2 *, for t>0, 

(5.3) 

(A -B)(A-T) = 2C^\ for t < , 

where we allow for the possibility of having a change in thermal entropy across the brane at 
t = 0, giving rise to different conserved quantities C + and C_ at the two sides. Continuity 
of the thermal entropy, C + = C- = C > 0, imposes B + V = 0. In this case the weak 
coupling constraint \B\ < A is automatically satisfied. The restriction of not exceeding the 
maximal temperature at the critical point, which is equivalent to A(0) = 0, further imposes 
B = T = and we obtain the non-singular, time-reversal invariant solution of section |4| 

We proceed to explore other solutions with C + ^ C_, preserving, however, the maximal 
temperature constraint, which requires B — Y. A very intriguing non-singular solution arises 
when C_ = 0, which also saturates the weak coupling constraint, A = B = T: 

p 2<t>o P A[\t\+t] 

e 2 <t> = e 2A = e 2A ° - (5 4) 

l + A[\t\+t] ' i + A [\t\+t)- l ° } 

For t < 0, the spacetime metric is flat, A = const., and the dilaton is constant. For t > 0, 
the Universe is filled with thermal radiation. The physical relevance of this solution can be 
understood as follows. The flat t < region, including the brane at t = 0, admits a Euclidean 
interpretation. The Euclidean path integral over it can be used to define the wave-function 
for the cosmology at t > 0, which summarizes the very early-time history of the Universe, 
while for later times the Euclidean brane creates the radiation. 

To understand this, we consider the Euclidean solutions, whose integration constants 



satisfy Eq.(5.3) after flipping the signs of the right hand sides. For C_ = 0, C + > 0, the 
solution maintaining A = B satisfies 2A(A + V) = -2C + e 2<t>0 = -(2A) 2 . The last equality 
follows by identifying the brane tensions in the Euclidean and Lorentizian solutions and 
yields T = —3 A. In half of this solution, for Euclidean time tE = z < 0, the dilaton is 
constant and the metric takes the form: 

ds 2 = e 2X °- 4Az [dz 2 + dx 2 ] = dp 2 + p 2 dy 2 where 2Ap = e x "~ Az , y = 2Ax . (5.5) 

In terms of the polar-like coordinates p and y, the metric is locally flat, with the radial 
coordinate p being bounded from below: p > e x °/(2A). When y — > iu>, the metric takes 
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the static Rindler form, where u plays the role of time. Although this space is locally flat, 
it is also intrinsically thermal possessing an effective "geometrical" (entanglement) entropy 
Thus, even though C_ is zero in this system, we can associate to the half-solution a non-trivial 
entropy The Rindler geometrical entropy is transformed via the brane into the thermal 
entropy of the late-time cosmological evolution. The half-Euclidean solution z G (—00, 0], 
which contains the brane, can be used to define the wave-function of the Universe described 



by Eq.(5.4) at t — 0. This mechanism will be investigated further in future work. 

Another interesting Euclidean solution, with C + = C_ = 0, corresponds to A = B = —T, 
including the brane localized at Euclidean time = z = 0: 

20o A[\z\-z] 
P 20 _ « 2A _ 2A « 

l+A[\z\+z] ' l + A[\z\+z}- 1 j 

As before, when £g = z < 0, the dilaton is constant and the metric is flat, taking the 



form given in Eq.(5.5) under the replacement A — >■ A/2. On the other side of the brane, 
tE = z > 0, both the metric and the dilaton become non-trivial 

p20 O a 7 i I 

e 2 * = - , ds 2 = e 2X ° aZ + f . (5.7) 

1 + 2Az l + 2Az y ' 

After an appropriate coordinate change, this solution becomes: 



e 



24> 



e 20o+2A o l l e 2A 

ds 2 = dp 2 + — dy 2 where Ap = e x °\fi + 2Az , y = —— x . (5.8) 



A 2 p 2 ' r p 2 * r ' * A 

In this parametrization, we recognize the well-known T-dual geometry of the two-dimensional 
flat space. Here, again, the coordinate p is bounded from below, p > e x °/A, thus avoiding the 
well-known singularity of the dual-to-flat metric at p = 0. The string coupling is bounded 
from above by its value at the position of the brane, g 2 tr < e 2 ^ . The dual-to- Rindler 
space is obtained by y — > iu, and so in this sense, it is also thermal possessing a non-zero 
"geometrical" entropy. On either side of the brane, the Euclidean solution is supersymmetric 
[32] , as is also implied by the choice C± = 0. The brane, however, breaks supersymmetry 
locally and creates the geometrical entropy in the Euclidean. The Euclidean solution clarifies 
further the importance of the brane. Namely, the brane glues the two dual geometries 
together in a consistent way. 



The time-reversed solution of Eq.(5.4) has C + = and C_ = C; it can be obtained when 
B = —A < 0. In this case the spacetime is flat and the dilaton is constant for positive 
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times t > 0. The Universe is filled with thermal radiation at negative times. As a result, the 
Universe is static with a thermal Rindler interpretation at later times. In this cosmological 
solution, the thermal entropy at early times is transfomed via the brane at t = to the 
geometrical entropy associated to the late-times Rindler space. 

Another interesting solution is obtained by combining these two solutions, assuming 
branes localized at two different times, namely at t = ±a. A time- reversal invariant solution 
can be easily derived, 

20o „A[|t+a|+|t-a|-2a] 
e 2</> = Z e 2A = e 2A « /g g\ 

l + A[\t + a\ + \t-a\-2a] ' 1 + A [\t + a\ + \t - a\ - 2a] ' y ' J 

with the following behavior: 

i) for t < —a, the solution for the dilaton and the metric is non-trivial, describing a con- 
tracting thermal Universe, with entropy coefficient Ce 2< ^° = 2A 2 . 

ii) for t > a, the dilaton and the metric have non-trivial time-dependence, describing an 
expanding thermal Universe with the same entropy coefficient as in i), Ce 2 ^ = 2A 2 . 

iii) in the intermediate region —a < t < a, the metric looks locally flat and the dilaton is 
constant. 

The intermediate flat region t G [—a, a], including the branes, admits a Euclidean inter- 
pretation, and as we explained before in the half-solutions it possesses a geometrical entropy 
and temperature in the Rindler cosmological frame. Thus the thermal entropy of region i) is 
transformed via the first brane at t = — a into the geometrical entropy of the intermediate 
region, which in turn is transformed back to thermal entropy in region ii) via the second 
brane at t — a. In the limit 2a — > 0, we recover the solution presented in section|4j Thus, the 
parameter a, describing the separation of the branes, provides a new relaxation time scale 
for the contracting-to-expanding cosmological phase transition. The whole picture can be 
understood via Euclidean instanton transitions between contracting and expanding thermal 
Universes. The pure Euclidean as well as the Rindler and dual-to-Rindler solutions will be 
investigated further in future work. 

6 Conclusions 

In this work we identified string-theoretic ingredients which can be used to resolve both the 
Hagedorn instabilities of thermal string theory, as well as the initial "Big Bang" singularity 
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of the induced cosmological evolution. In order to incorporate these ingredients, it was nec- 
essary to work beyond the low energy effective field theory approximation, examining purely 
stringy phenomena. For technical reasons, we focused on simple (4, 0)-supersymmetric string 
models, where adequate quantitative control over the dynamics is available. More specifi- 
cally, we considered the class of two-dimensional Hybrid models, which are characterized by 
a very special asymmetry between the left- and right-moving sectors. As in the usual super- 
string models, the left-moving sector is supersymmetric, giving rise to the 16 real spacetime 
supercharges. The right-moving sector is non-supersymmetric at the massless level, but it 
enjoys the MSDS symmetry structure, which guarantees boson/fermion degeneracy in all of 
the right-moving massive levels. 

We considered tachyon-free, thermal configurations of the Hybrid models, where in ad- 
dition to the temperature, we turn on certain discrete "gravito-magnetic" fluxes threading 
the compact Euclidean time circle. These fluxes inject non-trivial momentum and winding 
charges into the thermal vacuum, which effectively regulate the contributions to the thermal 
partition function of the massive string oscillator states and restore the thermal T-duality 
symmetry of the finite temperature stringy system. As we have seen, the tachyon-free ther- 
mal configurations can be described in terms of freely acting asymmetric orbifolds, which act 
thermally on the left-movers but leave the right-moving MSDS symmetry unbroken. Thanks 
to the latter symmetry, the one-loop thermal partition function was explicitly calculated, be- 
yond any a '-approximation. 

The thermal Hybrid systems have three characteristic regimes, associated to the space 
of light thermal momenta, the space of light thermal windings and the extended symmetry 
point, which arises at the self-dual value of the thermal modulus and connects the two. 
At the extended symmetry point, extra massless thermal states appear, with a clear brane 
interpretation in the Euclidean. These localized states are responsible for a phase transition 
between the two dual thermal phases. The physical temperature of the system is bounded 
from above, attaining its maximal critical value at the phase transition. 

An important result of this work is that we succeeded in writing down a stringy effective 
Lorentzian action able to incorporate the three stringy regimes simultaneously, which is 
exact in a 1 at the two derivatives level and is valid up to genus- 1. The essential feature 
in this action is the spacelike brane contribution that glues together the spaces of light 
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thermal windings and light thermal momenta. This brane is sourced by the additional 
massless thermal states, which appear at the extended symmetry point. We explicitly showed 
how the sign and the value of the brane tension are determined in terms of the "rapidity" 
coefficients characterizing the classical backgrounds of the corresponding localized thermal 
fields. The spacelike brane, together with the bulk thermal corrections, induce a non-singular 
cosmological evolution, describing a bouncing thermal Universe. The bounce coincides with 
the phase transition between the two dual Hybrid thermal phases. The cosmology remains in 
the weak coupling and low-curvature regimes, provided that the rapidity coefficients which 
determine the tension of the brane are small. Continuity of the thermal entropy across the 
brane uniquely determines the cosmological solution, which is time-reversal invariant. 

Separated spacelike branes, localized at different times, transform the thermal entropy 
into a geometrical entropy associated with a quasi-static thermal Rindler space (or its dual) 
and vise-versa. Utilizing this property of the branes, we obtained more general, time-reversal 
invariant solutions, which are also free of the initial Big Bang singularity. The obtained 
non-singular solutions describe transitions between a contracting thermal Universe and an 
expanding one. The effect of the spacelike branes is to induce an intermediate relaxation-time 
period during which the scale factor of the Universe is quasi-static, and the entropy appears 
in a geometrical Rindler-like form. The intermediate locally flat region, including the branes, 
admits a Euclidean interpretation describing instanton transitions among contracting and 
expanding thermal Universes. We would like to investigate this connection further in the 
near future. 

To our knowledge, the two-dimensional thermal Hybrid model presented in this work is 
the first example in the literature where both the stringy Hagedorn singularities as well as 
the classical Big Bang singularity are both successfully resolved. 

We believe that phenomenologically relevant, higher dimensional string models exist, 
where the essential ingredients, used in the thermal Hybrid model to overcome the initial 
singularity problem, are also present. Namely: i) the restoration of thermal duality symme- 
try in the presence of non-trivial discrete "gravito-magnetic" fluxes, ii) the appearance of 
new light thermal states at the self-dual point, giving rise to a localized higher- dimensional 
brane interpretation and iii) the gluing of dual thermal phases via higher dimensional space- 
like branes, before the occurrence of the Big Bang singularity that is always present in the 
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naive field theory approximation. We conjecture that in the higher dimensional cases, the 
resulting duality structure which encompasses all these ingredients will be much more in- 
tricate, involving other supersymmetry breaking scales and perhaps non-perturbative string 
phenomena as well. Although the technical challenges in the higher dimensional cases are 
expected to be more delicate, the key property, which is the joint effect of string dualities 
and branes "protecting" the stringy system from the occurrence of spacelike cosmological 
singularities, could still be realized. Obviously, more work needs to be dedicated to this 
ambitious direction in order to obtain non-singular cosmological models that will be phe- 
nomenologically viable. 
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